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ABSTRACT 

Gamma-ray bursts have the potential to produce the particle energies (up to 10^^ eV) 
and the energy budget (10^^ ergyr"^ Mpc~'^) to accommodate the spectrum of the high- 
est energy cosmic rays; on the other hand, there is no observational evidence that they 
accelerate hadrons. The Fermi GST recently observed two bursts that exhibit a power- 
law high-energy extension of the typical (Band) photon spectrum that extends to ~ 30 
GeV. On the basis of fireball phenomenology we argue that they, along with GRB941017 
observed by EGRET in 1994, show indirect evidence for considerable baryon loading. 
Since the detection of neutrinos is the only unambiguous way to establish that GRBs 
accelerate protons, we use two methods to estimate the neutrino flux produced when 
they interact with fireball photons to produce charged pions and neutrinos. While the 
number of events expected from the Fermi bursts detected to date is small, we conclude 
that an event like GRB941017 will be detected by the IceCube neutrino telescope if 
gamma-ray bursts are indeed the sources of the observed cosmic rays. 

Subject headings: gamma rays: bursts — gamma rays: observations — neutrinos 



1. Introduction 



The sources of the extragalactic cosmic rays with energies in excess of ~3xl0 eV remain a 
mystery, but one of the best motivated candidates is gamma-ray bursts (GRBs). Large cosmic-ray 
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energies can be achieved in the prompt phase of the G RB firebah where internal s hocks have the 
potential to accelerate charged particles up to ~10^^ eV (jVietrilll995l . IWaxmanll 19951 ). Additionally, 
the total energy density in the Universe of cosmic rays must be matched by a sufficiently high 
hadronic energy density in the GRB with Phad.GRB ~ Pen- GRB observations identify synchrotron 
photons produced by the electrons accelerated in the fireball with an energy 6gE!^Qrp ^ 10 ergs, 



TOT 

where EISj^Q^p is the total isotropic energy released by the burst. GRB fireballs also carry energy 
e^EJI^Q^ in the form of magnetic fields, and, if they are the sources of cosmic rays, energy e^EJ^Q^ in 
non-thermal protons. Further energy is carried by the thermal leptonic and hadronic parts, which 
is not discussed here. 

GRBs emerge as credible sources for the ultra high-energy cosmic rays because their observed 
flux can be accommodated with an energy density in protons that is similar to that in electrons. 
Recent estimates of the local rate of GRBs yield a maximum of no~l Gpc~^ yr"^ assuming that 
GRBs follow the star forma tion rate. For a st ronger evolution with redshift, the local rate can be 
as low as 0.05 Gpc~'^ yr"^ ( Asano et al.ll2009l ). Using this result we estimate the electromagnetic 
energy density from GRBs to be in the range /0<,ni,GRB ~ ho eg ^tqt = 5 x 10^^ - 10^^ ergs Mpc~^ yr~^- 
In order for GRBs to be the sources of cosmic rays, their hadronic energy density p^ad grb needs to 



10 ergs Mpc 



produce the observed cosmic-ray energy density Phad. grb 
theref ore c onclude that e„/ee ^ 1-20. Discussion in, for example. IWaxman et al\ (|l997l ). 



yr 



We 



Waxman 



(j2004l ). and lMurase et al\ (j2008l ). agrees that the ratio must lie around those values. As it can be seen 
from the rough estimate above, the actual number strongly depends on the local rate of GRBs, 
which is still quite uncertain. In addition, the non-thermal extragalactic spectrum is expected 
to extend to energies below the knee, but reni ains unobserved due to the larger contribution of 
galactic cosmic rays (see e.g. lAhlers et al\ (|2005l )). If this is the case, the average fraction of proton 
to electron energy needs to be larger. 

The Fermi gamma-ray satellite observatory recently observed two bursts, GRB090510 and 
GRB090902b, that show a st atistically significa.nt deviation from the typical GRB spectrum de- 
scribed by the Band function (jAbdo et al\ l2009l : iBissaldi et al\ l2009l ) . A flux of high energy events 
is detected that extends to energies of ~30 GeV following a power-law spectrum (T able [J). This is 



simil ar to the much more luminous burst observed by EGRET in 1994, GRB941017 (jGonzalez et al. 



20031 ). There has been much discussion on the possible origin of these high-energy tails. A lep- 



(2003): 


Pe'er et al. 


J2OO4 


); 


Stern et al. 


(2004) 



also connected to the emission of high energy photons. Prot on synchrotron emission is discussed in 



(jPermer et al.ll2004l : iRazzaque et adl2009l : lAsano et a/.ll2009l ). while the production of neutral pions 
in photohadronic reactions , which lead to the emission of h igh-energy photon s, is discussed in e.g. 



Alvarez-Muhiz et al\ (|2004l ) for GRB941017. It is noted by lAbdo et al\ (|2009l ) that at least in the 
case of GRB090510, all models show difficulties in explaining the emission and there is as yet no 
model that perfectly fits the observation. In this paper, we will investigate the possibility that the 
high energy component results from vr'^-decays and is therefore a signature of proton acceleration. 
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In particular, this model requires relatively high baryonic loading of the jet. We will discuss this 
fact in detail in this paper. 

The main contribution to the initial opacity of the fireball comes from the annihilation of 
photons into pairs. The Fermi observation of a non-thermal spectrum up to an energy i?max of 
tens of GeV can be used to constrain the minimum bulk Lorentz factor Tmin required to make the 
source optically thin at the time of the gamma-ray display. For all photons with energy E < -Bmax 
the condition t^^{E) < 1 must be fulfilled where r^^ is the opacity. The observation of photons 
with energies of tens of GeV requires highly relativistic outflows with F ~ 10^. Because, on the 
other hand, the observed energy fiux of order 10"'^ ergscm^^ s^^ is typical of an average burst, the 
large boost factor implies that the photon density in the rest frame of the burst is low. This is a 
strong effect as the photon density is suppressed by F~^. From t^^{E) = 1, we can determine eg 
by finding the electromagnetic energy over the volume of the fireball as a fraction of the total GRB 
energy. The optical depth is defined as 

^11 = —- (1) 

Here, AR is the thickness of the fireball shell in its rest frame and A^^ is the photon mean free 
p ath. From the definition of the mean free path and the photon number density as given e.g. 



in iGuetta et al\ ( 20041 ) . we then obtain 

T^^ = ARa^^ = ARa^^ [ ) = ARa^^ ( ,^_^,7^\.^^ ) • (2) 



Therefore, since the isotropic luminosity = eei?xoT/^90i the condition Try^{E) = 1 gives 

-2 r \ V 5t \ ( E^ \ f Tgo\ f 10^3 ergs 



where F is the bulk Lorentz factor of the fireball, St is the variability timescale, E^ is the charac- 
teristic gamma-ray energy (which we take to be the peak energy of the event), Tgo is the duration 
of the burst, and -E^ot is the total isotropic energy of the burst. We therefore estimate Ce for 
GRB090510 to be ~ 0.05 and for GRB090902b to be ~ 0.02. From the low values of ee thus 
obtained protons seem to dominate the fireball. However, these numbers depend on having input 
a value for -Etot' which cannot be measured, and so precise values of eg cannot be determined. 

In this paper we will first discuss the properties of the bursts. We subsequently compute the 
neutrino fiux inevitably produced when the protons interact with fireball photons. Their obser- 
vation would provide incontrovertible evidence for the pionic origin of the additional high energy 
component in the burst and support the speculation that GRBs are the sources of the highest 
energy cosmic rays. 



Is a kilometer-scale neutrino telescope such as IceCube sensitive enough to shed light on these 
questions? High energy neutrinos are produced in the fireball when protons produce pions in 
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interactions with the photon field. Using the A-resonance approximation, 



PI 



A+ 







1/3 of the cases 
2/3 of the cases 



(4) 



This gives pion ratios of vr^ 
pions decay as tt"*" — > ^"^ v 



TT^ = 1 : 2. The neutral pions decay as vr'^ 27, and the charged 

-T' I'ei^fj.i'fi- Here, a single neutrino carries approximately 1/4 of 




the TT"*" energy and a photon carries 1/2 of the vr energy. The calculation of the neutrino flux 



Becker et al 



200J with the following results: whereas an average burst produces only ~ 10 



(IWaxman et al\ 119971 1 has been performed in detail for the BATSE bursts (jGuetta et al 



2004 



neutrinos, bursts that are unusually energetic or nearby may produce an observable flux in a 
kilometer-scale neutrino telescope of order 10 events per year. We suggest that the power-law 
high-energy spectral feature can identify such bursts. 

We will compute the neutrino fluxes expected in IceCube using two methods: the standard 
fireball model, and the bolometric method which relates the energy in neutrinos from the decay of 
charged pions to the observed photon energy assuming that it is of pionic origin. We will conclude 
that while the neutrino rates from the Fermi GST bursts are unexceptional, a burst like GRB941017 
extending to tens of GeV energy will be detected by IceCube. IceCube observes cosmic neutrinos 
in a background of neutrinos produced in the atmosphere. Given that neutrinos of GRB origin are 
relatively energetic and that the direction and time of the events can be correlated to satellite alerts, 
the atmospheric background is suppressed and few neutrino events may still represent a conclusive 
detection. We note that relating neutrinos to an observed gamma-ray flux using the A-resonance 
approximation is extremely conservative in thi s context and the po tential neutrino flux could be as 
much as a factor 4 greater (jRachen et aLlll998l : iMurase et aLll2006l ). 



2. Detection of GRBs at high photon energies 



The main contribution to the opacity of GRBs fireballs comes from the annihilation of pairs 
of photons into pairs. The observation of a non-thermal spectrum with maximum energy E'max 
requires the fireball to be optically thin to photons at the ti me of emissioii and t herefore can be 
used to constrain the bulk Lorentz factor T. By generalizing iLithwick et al\ (j200ll ) to account for 
arbitrary photon spectra and the energy dependence of the pair production cross-section, we can 
express the optical depth r for the most energetic photon in the GRB fireball (with observed energy 



as 



AndlSt 
47r(r2c<5t)2 



d(cos0) 



(1 - cos9) 



f 

Jo 



(l + z)E^ (l + z)En 



COS ( 



(5) 

where dL{z) is the lumin o sity d istance to the source assuming the same ACDM cosmological param- 
eters used in lAbdo et al\ (|2009l ). a^-y{Ei, E2, cos 9) is the cross-section for tw o real photons col liding 
at an angle 9 with energies Ei and E2 to produce an electron-positron pair (jCould et a/.lll967l ). and 
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dN^/dE^ is the observed photon spectrum (on Earth) from the GRB. The colhsion kinematics must 
take into account the boost factor of the fireball T and the fact that the observed photons have 
nonzero redshift z. The additional factors of F and (1 + z) modifying the arguments of a^^ there- 
fore transform the observed photon energies to the energies in the center-of-mass frame and at the 
source, respectively. We consider only the contribution to the opacity due to pair production and 
neglect other processes such as inverse Compton scattering and synchrotron self -absorption. The 
opacity due to Compton scattering is small in a GRB fireball (iRazzaque et aLll2004l) and sy nchrotron 
self-absorption is not important at photon energies below ~ 10^6 eV (jMurase et al\mm . 



The first clear detection of a GR B with a power-law component in addition to the standard 
Band spectral form (jBand et al\ Il993l ) was made by the EGRET satellite with the detection of 
GR B941017 with a spect rum that extended to 200 MeV, reaching the hmit of EGRET's sensitiv- 
ity (jGonzalez et aLll2003l ). Estimates of the correspon ding neutrino flux and the expected event 
rate in km^-scale neutrino telescopes are presented in lAlvarez-Mufiiz et al\ (j2004l ) . Of particular 
importance is the fact that since EGRET did not detect a cutoff of the power-law component, the 
flux can be modeled as extending to potentially very high energy i^max- Using Eq. [5] with the total 
observed photon spectrum (Band + power law), we may then find the boost factor T for which 
T = 1. We call this Tmin, the minimum boost factor for which the fireball is transparent to all 
photons up to energy -Bmax) and is shown in Figure [1] for time bins 2-5 (time bin 1 showed no 
significant power-law component). 

The Fermi Gamma- Ray Space Telescope has been collecting very energetic photons (8keV- 
300 GeV) from gamma-ray bursts since June 2008. As of September 2009, at the end of the first 
year of observations, 10 bursts have been detected by the GBM (Gamma-ray Burst Monitor) and 
the LAT (Large Area Telescope) simultaneously with photon energies above ~100 MeV. Fermi LAT 
detections have constrained the values of the jet boost factors to values of F ~10^ for GRB080916c, 
GRB090510 and GRB090902b thus confirming the highly relativistic nature of the outflows. Fur- 
thermore the time integrated spectra of the recently observed GRB090510 and GRB090902b show 
statistically significant deviations from the Band function that extend up to ~30 GeV. 



GRB090510 is a short burst (Tgo 
up to ~30 GeV with a fluence of 



2.1s) that shows a spectral high-energy component 
(~ 40% of the total fluence in the en- 



~1.8xl0 ^ ergs cm ^ 



ergy range 10 keV- 30 GeV) (jAbdo et al\ l2009l ). The isotropic energy release is estimated to be 



ee-Eiso ~1.08x 10^'^ ergs. The ~31 GeV photon is detected ~ls after the trigger time and sets 
the highest lower limit on a GRB Lorentz factor (F ~1200) thus proving the outflows in short 
gamma-ray bursts to be as highly relativistic as those in long gamma-ray bursts. 

The extra component of the long (Tgo =21.9 s) GRB090902B extends up to 11.2 GeV and 
its fluence accounts for ~24% of the to tal fluence over the energy range lOkeV-lOGeV in the 
first 25 seconds of the prompt emission (iBissaldi et aLll2009l ). The corresponding isotropic energy 
release is measured to be eg -Eiso ~3.63x 10^^ ergs. The delay of the highest energy photon with 
respect to trigger time is ~80s and the Lorentz factor of the jet is estimated to be F ~1000 
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Fig. 1. — For time bins 2-5 of GRB941017, values of F such that r^^ = 1 for a photon of energy 

-^max (-Tmin)' 



from opacity considerations involving the highest energy observed photon (11.16 GeV) during the 
prompt emission phase. In Fig. [21 a comparison between the fluences of the Fermi bursts and the 
fluence of GRB941017 as a function of i?max is presented and the relative weakness of the Fermi 
bursts appears evident. In Table [H the spectral parameters from the Fermi GST observations of 
GRB090510 and GRB090902b are summarized. 

Deviations from the Band-only fit in the spectra are particularly interesting in the context 
of hadronic acceleration within the fireball and relate closely to predictions of neutrino fluxes de- 
tectable on Earth with km^ telescopes. If these extra-components originate from vr'^- decay photons 
they provide an optimal benchmark for testing models of hadronic acceleration in GRB engines. 



3. Fireball neutrinos 



In the hadronic fireball a burst of high-energy neutrinos is expected to accompany the observed 
prompt fiux of gamma-ray photons. Assuming that electrons and protons are shock-accelerated in 
the same region, the neutrino spectrum can be calculated fro m the observed spe ctru m in gamma- 
rays using conventional fireball phenomenology as described in lCuetta et al. (2004) and Abbasi et al. 
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Fig. 2. — Fluence of the high energy component above 30keV as a function of the maximum energy 
£;max for GRB941017 and observed fluences for GRB090510 and GRB090902b. The arrow above 
the data point for GRB941017 emphasizes that EGRET was sensitive only to the part of the fluence 
below 200 MeV and that therefore the data point is a minimum. 



([20091). For a typical GRB, the gamma-ray spectrum is usually well described by a Band function: 



A- < 



100 keV I 

\ rv — R 

[a - (3)e-, \ 
100 keV 



E. 



exp(/3 — a) 



E^ 
100 keV 



if (q - /3)e^ ^ E^ 



if (a — I3)e-y ^ E^ 



(6) 



The interaction of accelerated protons with a power-law distribution with GRB photons results in 
a broken power-law neutrino spectrum. The resulting neutrino spectrum is then expected to follow 



dE^ 



'I 



e; 



if < ei 

if ei < E^ < £2 



(7) 



where the spectral indices are related to the indices of the photon spectrum according to 

Ou = 3- (3^ l3i^ = 3- a.y = l3i^ + 2. (8) 
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The first break energy ei is determined by the production threshold for the A-resonance, where 

. 1 / r \^ /MeV\ 

is the break energy of the Band function, T is the jet boost factor and z the redshift of the 
source. The spectrum steepens when pions lose energy due to synchrotron radiation prior to decay; 
the second break energy 62 is given by 



,7n^^r 1 f ^ V f St \ /I052ergss- 



l + z y €B yW^-^ J VlOms/V L 



^2 = W'GeY-—J- -— J . (10) 



7 



ee is the fraction of the total burst energy available to electrons, is the fraction of total energy 
going into the magnetic field B and 6t is the variability timescale. is the isotropic luminosity of 
the Band function, given by 

E- 4ndl{z)F- 



with At = Tgo and the luminosity distance for = 0.27, Q\ = 0.73 and h = 0.71. The Band 
gamma-ray energy fluence is defined as 

r°° dN 

F- = At dE,E,^ (12) 



7 



where the gamma-ray spectrum dN^/dE^ is the Band spectrum given in Eq. O 



Table 1: Spectral parameters of Fermi GRBs with power-law components. The total fluence in 
gamma rays F™^ is the sum of the Band fluence F]^ and the power-law fluence F™. The flu- 
ence for GRB900510 is calculated over the energy range 10 keV-30 GeV, and the fluence for 
GRB09092b is calculated over th e range 10 keV-10 GeV. The parameters for GRB941017 are 
tabulated in iGonzalez et al\ (j2003l ) . 





GRB090510 


GRB090902b 


z 


0.903 


1.822 




2.1s 


21.9s 


pTOT 

7 


5.02 X 10~^ ergs cm"^ 


4.36 X 10^^ ergs cm^^ 


7 


1.84 X 10^*^ ergs cm^^ 


1.05 X 10^^ ergs cm^^ 




-0.58 


-0.61 


/57 


-2.83 


-3.80 


^7 


2771 keV 


522 keV 


r. 


-1.62 


-1.93 




30.53 GeV 


33.40 GeV 


r 


1260 


1000 
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The normalization of the neutrino flux is determined by the efficiency of pion production. The 
relation between the gamma-ray and neutrino spectra is given by 



where 



1 €p 

X = - — 

See 



1-fl 



(13) 



and 



In Eq. m 



AR 



1052 



erg s" 



10ms\ / 102-5 
"IT 



(14) 



^tt) = 0.2 is the fraction of proton energy transferred to a pion in a single interac- 
tion and e^i is the break energy of the Band function. The neutrino energy fluence Fy is defined 
analogously to F^ (Eg. [T2|) . 



Table 2: Parameters for fireball neutrino fluxes. Typical values of the WB model are reported. 
Quantities marked with ' * ' are not determined for the particular burst and standard WB model 
values are used. 





WB 


GRB090510 


GRB090902b 




0.1-1 MeV 


2.771 McV 


0.522 MeV 




-1 


-0.58 


-0.61 




-2 


-2.83 


-3.8 


r 


300 


1260 


1000 




-1 


-0.17 


0.8 




-2 


-2.42 


-2.39 




-4 


-4.42 


-4.39 




10^-10^ GeV 


1.2 X 10^ GeV 


1.3 X 10^ GeV 


ev2 


10^ GeV 


3.6 X 10* GeV 


5.6 X 10* GeV 


z 


1-2 


0.903 


1.822 


it 


0.01 s 


0.01 s* 


0.053 s 




2-1000 s 


2.1 s 


21.9 s 



The values used to calculate the expected fireball neutrino rates from the FGST bursts, using 
only the Band component of the photon spectra, are listed in Table [2] together with the typical 
Waxman-Bahcall (WB) parameters. The numbers of neutrino events in IceCube estimated using 
the fireball formalism with the Band photon spectrum are reported in Table [3l 
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4. Neutrinos from power-law GRB spectra 

4.1. Fluences of secondary particles 

If we assume that the observed power-law components of the bursts discussed here are due to 
the decay of neutral pions produced in interactions between protons accelerated by shocks in the 
jet and photons, we can predict the accompanying flux of muon neutrinos from the production and 
decay of charged pions. The observed power-law spectral components are quite flat {E~^ — E~^'^), 
a spectral behavior that is compatible with synchrotron radiation from highly relativistic electrons 
and positrons which are produced when high-energy photons scatter in the fireball photon field. 
For interactions of protons accelerated in shocks with fireball photons, we assume that the final 
neutrino spectrum will be of the form of Eq. 7. If we find the charged pion spectrum that, upon 
decay, produces the correct neutrino spectrum, we can determine the corresponding neutral pion 
spectrum and thus the gamma-ray spectrum. Assuming that the cascading process conserves 
energy, we can normalize the fiuence of gamma-rays to the fluence of the measured power-law tail 
assuming its extension to some maximum energy -Emax 



E, dE, = E, dE„ (15) 



where £^min is the minimum measured photon energy and -Emax is the proposed upper limit of the 
measured power-law photon spectrum dN^^ /dE^. The upper limit of 10^^ eV for the uncascaded 
gamma-ray spectrum assumes that the parent protons extend to ~ 10^'^ eV and that gamma-rays 
take on average 1/10 of the parent proton energy. Due to the flatness of the measured power-law 
spectra the precise value of -Emin is unimportant and the energy going into hadronic gamma rays is 
determined purely by -Emax- Since we assume that the flreball is transparent for photons of energy 
-Emax; the boost factor F will vary as we vary the total energy in hadronic photons. 

The algorithm for flnding the flux of neutrinos is as follows. For a given -Emax; the measured 
parameters of the GRB will determine the minimum boost factor Fmin (Eq. 5). These parameters 
and boost factor determine the neutrino spectral indices and break energies (Eqs. 8-10), giving us 



Table 3: Fireball neutrinos in IceCube. The number reported here for GRB941017 is the sum of 
the number of events from e ach individual t i me b in. The parameters of the photon spectra for 
the time bins are taken from I Gonzalez et al\ (j2003l ) and parameters from Table [2] are used in the 
calculation of the number of neutrino events. 





number of fireball z^'s 


GRB941017 


9.7 X 10^2 


GRB090510 


1.6 X 10""^ 


GRB090902b 


1.9 X IQ-'^ 



-li- 



the unnormalized neutrino spectrum. Assuming that the neutrinos take on average 0.25 of the pion 
energy, we define the charged pion spectrum to be 



( el';-^^ E-^'^ if^. <ei,. 

E^^" if 61,. <E^< €2,^ (16) 

el^;^^E-'<^ if E^ > £2,. 



with pion break energies ei^-j^ relating to the neutrino break energies (Eqs. [Qland llOp via 64,,^ = 4ei. 
From the charged pion spectrum, we derive the neutral pion spectrum. 



dE^ 



dN^+ 

II £/7r < e2,7r 



We assume that the photohadronic interaction takes place at the A-resonance, giving twice as many 
neutral as charged pions. Moreover, the second break in the neutrino and charged pion spectra is 
due to the cooling of long-lived charged pions in the fireball and will not be present in the neutral 
pion spectrum. 

From the pion spectra we can find the decay spectra of gamma rays and neutrinos, normalized 
relative to the arbitrary factor of /^+ in Eq. [16] (see Fig. [3|). Integrating over the neutrino and 
gamma-ray spectra then gives the amount of energy going into neutrinos relative to gamma rays. 

K^dE^ = rj E,^dE„ (18) 

CL^u Jo "-f^7 

The procedure for obtaining the decay particle spectra from the pion spectra Eqs. [16] and[l7] is 
described in Section 4.2 below. Finally, using Eq. [15] with the factor r] found from Eq.[T8l we can find 
the actual amount of energy in the gamma-ray spectrum and therefore the absolute normalization 
of the neutrino flux: 

Since neutrinos originating from pion decay have the flavor ratio v^. : : Ur = ^ '■ '■ over an 
astronomical baseline oscillations will transform the beam into the flavor ratio 1:1:1 if ^13 is small. 
Therefore the flux of muon-type neutrinos that reach the detector is half of the emitted flux and 
the final neutrino spectrum must be multiplied by a factor 0.5. 

In this method it is clear that higher observed maximum gamma-ray energies lead to more 
energy in neutrinos. This does not contradict the discussion in Section 2 that describes how a 
higher -Emax implies a larger boost factor F and hence a more transparent fireball. While a larger 
boost factor will imply a smaller photomeson production efficiency, relating -Emax to the proton 
fraction ep/ee via Figure [I] and Eq. [3] gives us Cp/ee ~ 10,000 for -Emax ~ 50 GeV, to keep the 
total energy of the burst constant. A proton loading of this magnitude would likely compensate for 
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Fig. 3. — Pion, neutrino and uncascaded gamma-ray spectra for Bin 2 of GRB941017, assuming 
a boost factor T = 300. The spectra are not normalized but do account for the effect of neutrino 
oscillations. 



the lowered photomeson production due to increased boost factor. The question of the baryonic 
loading is discussed in detail in Section [5j 

We also assume in this section that the observed power-law components correspond to all the 
energy emitted in hadronic gamma rays. There is, however, the possibility that some part of the 
high-energy gamma-ray tail will escape the GRB without interaction if the opac ity is not large 



These photons would ca scade outside the source and are typically not observed (jRazzaque et al. 



20041 : iMurase et a/.ll2008l ). In that case, estimating the total hadronic gamma-ray energy using the 
observed photons would result in an underestimation of the energy and hence the neutrino flux. As 
a result, the neutrino fluxes derived in this section are from a minimum hadronic energy and are 
therefore conservative. If GRBs are optically thin to photohadrq nic interactions, synchrotron self 



Compton (IGranot et al. 



2003 



Pe'er et al. 2004; Stern et al 



20041 ) or proton synchrotron radiation 



(jPermer et al\ l2004l : iRazzaque et al\ l2009l : lAsano et al\ l2009l ) could dominate. Here, however, we 
posit that the sources are not optically thin to photohadronic interactions. It is clear that the 
observation of neutrinos with a full IceCube detector will resolve this issue. 



- 13 - 



4.2. Pion decay spectra 



The muon neutrino and gamma-ray spectra are determined from pion spectra as follows: Both 
gamma rays are the product of the decay of the neutral pion tt" — >■ 7 + 7. The two muon neutrinos 
are due to the decay chain vr^ v ^ + /x^ , /x^ — )• ^ + Ve + e^. 

The Kamma, - rays and the first neutrino are found from standard two-body decay kinemat- 
ics (lDermerill986l : iGaisserlllQQOl ): 



min (-^z ) 

For both particle types {i = 1^,7), 



E, 



1 



(20) 



(21) 



where = (neutral pion decay), r^, = {m^/m-n-)'^ (charged pion decay), = 2, and A,^ = 
(1 — r^)~^. Due to the E + E~^ form of -Emin, the two-body spectra are symmetric (on a log-log 
scale) around a peak set by the pion mass and by the decay kinematics scale factor r. This peak 
is at m7r/2 ~ 70 MeV for neutral pion decay and (1 — rjy)m,r/2 — 30 MeV for charged pion decay. 
For a power-law distribution of pions with energies much larger than the pion mass the spectrum 
is a power-law of the same slope. 

The second neutrino is due to the decay of the muon from the charged pion fi^ ^ -|- Pe + i^ii- 
T his is a three-b ody decay of a particle with a two-body decay energy distribution. The spectrum 
is dCaisserllioOol ) 



''>U,2{E,y 



dE„ 



E-n I 



dE^ 



b^^lE^^) 1 dn 



i?Mmin JE^^,,, (1 - r^) ^El - ml E^ dy 



The full limits on the integrals, valid at all energies, are 

,2 



E, 



fi mm 



Eu + 



4.E^ 



Et, 



E-Kl 



'Ue,. 



(22) 



Gaisseij (|l990l ) considers only the spectrum at high energy and therefore omits the second term of 
each limit. The last term in the integrand {1/ E^)dn/ dy is the muon decay distribution. It is given, 
for y = E„/E^, by 



dn 1 

dy Jx 



2y 



dx 

X 



(23) 
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where 



/o = 2x2(3 -2x) 

/i = 2x^(1 -2x) 

1 / 2E^r^ 



ibP^ corresponds to the decay of /j,^ , respectively. The hmits of the integral are 

22/ 



1 + /3m 

2y 

min[l, 



1-/3. 



As a practical matter, these spectra involve integrals to +oo. This means that a change of 
variables is required before the integrals can be evaluated numerically. In this work we used the 
transform 

'^Z" „ /1\ dt 







f{x)dx= I f{-j^- (24) 



We note in particular that the more general transform 



fix)dx= r/fa + ^^^^ 







does not evaluate the integrals for pion decay correctly around the low-energy kinematic peak. In 
the case of muon decay (Eq. [22]) the integral transformed with Eq. [25] does not converge. 



4.3. Neutrino spectra and event rates 

Folding the derived neutrino fluxes with the IceCube effective area, we find the number of 
neutrino events as a function of neutrino energy and -Bmax of the observed gamma-ray fluence. 
In Fig. d] we show the number of detected neutrinos with energy greater than detector threshold 
as a function of i^max- We call attention to the fact that a burst with the same parameters as 
GRB941017, extending to i?max > 50 GeV, will produce > 1 neutrino event in IceCube. Given the 
absence of background events over such a short time interval, only a few events would be needed 
to constitute discovery of proton acceleration in GRBs. While very high E'max would require a 
(possibly) unphysically large boost factor and proton fraction, the detectability of a GRB941017- 
like burst with -Emax within Fermi's energy range is highly encouraging. 

These neutrino event numbers are very different from the numbers derived with the standard 
fireball phenomenology using only the Band spectrum (Table [3]). While the two numbers of neutrino 
events for GRB941017 are similar if -Emax ~ 200 MeV, the fireball event numbers for the Fermi 
bursts are both more than an order of magnitude smaller than the numbers calculated using the 
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Fig. 4. — Neutrino events in IceCube from 3 GRBs as a function of maximum observed gamma-ray 
energy. 

bolometric method, showing the large contribution of the GeV power-law component gamma rays 
relative to the lower-energy Band spectrum photons. 



5. Photomeson production efficiency 

The dependence of the pair-creation opacity r^^ on the jet boost factor T was presented in Eq.[2] 
where it was shown that r^^ oc F"^. The source must be therefore more transparent when photons 
with increasing energy E'max are detecte d. In addition to thi s, the efficiency of pion production. 



was shown to be proportional to r^^ (jWaxman et a/.l 119971 ). The predicted flux of neutrinos is 



therefore believed to be suppressed when F becomes large. Here, we argue that this problem can be 
overcome by the requirement that GRBs are proton dominated. The parameter Cp/ce (an intrinsic 
burst parameter, hence a constant during the evolution of the burst) has to be ^ 1 to compensate 
for the decrease of the efficiency of pion production at large boost factors. However, since this is 
only the high-energy, nonthermal part of the energy balance, extremely high values of Cp/ce cannot 
be ruled out. 
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In our discussion we have assumed that the high-energy component of all time bins of GRB941017 
originates from the interaction of UHE CRs with photons of the spectral region described with the 
Band function. Here, we sum up the fluences of time bins 2-5 to get the total (time-integrated) 
fluence of the observed high-energy component: 

5 

-^HE^(-^max) = y^^HE(-^max) (26) 
i=2 

where: 

-^HE(-^max) = / E^—^dE^ (27) 
JSOkoV "-f^7 

and dNi^^/dE^ is taken from Gonzalez et al.(2003). The index i runs over time bins. The value 
of -Emax is treated as a free parameter in all time bins since no break was detected by BATSE up 
to 200 MeV. Therefore, we let it vary up to 1 TeV and we calculate the corresponding observed 
total fluence of the high-energy power law. We then compare F™'^(ii^max), which is proportional 
to ep/ee- 

F™- ^ (28) 

with the theoretical flux of photons from vr'^-decay, -Fhe°"^- this way the hadronic model can 
reproduce the observed high-energy component for a given value of ep/ee- In Fig- Ethe ratio ep/ee 
is plotted against i?max- It is important to stress that the correct relation between T and -Emax 
has been used. We also note that for this GRB the values of ep/ee shown in Fig. [5] refer to the 
time integrated fluence defined in Eq. 26. This is due to the assumption that ep/ee is an intrinsic 
property of the burst and that therefore it should not vary with time. 

For each value of ep/ eg an estimate of the beamed energy release is then inferred using the relation: 

-E'bcam ~ ^^^2 

Results are shown in Fig. [6] and a physical limit can be inferred from the figure. Very recently pair 
instability SNe have been shown to release > 10^^ ergs (jA. Gal- Yam et al.ll2009l ). A similar energy 



release can be realistic in GRBs. Thus, we conclude that if beaming effects are considered the total 
energy release implied by the GRBs discussed here can be accomodated. 



Applying the same argument to the time integrated spectra for GRB090510 and GRB090902b up 
to the observed -©max yields ep/ee ~ 170 and ~63 respectively and the beamed energy release using 
Eq. [29] with the corresponding numerical value of the jet boost factor (table [2]) is estimated to be 
~ 3.66 X lO^s ergs and ~ 8.70 x lO''^ ergs (Fig. E]). 
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Fig. 5.— GRB941017: e^/ce Vs. observed E„ 



6. Summary and Conclusions 



Three gamma-ray bursts with statistically significant high-energy power-law spectral compo- 
nents have been detected thus far. While both leptonic and hadronic models have been proposed to 
explain the additional components, all models show some difficulties in explaining the observations. 
In this paper, we discuss the possibility of proton-dominated bursts with photohadronic interactions 
being responsible for the production of the high-energy component. We calculate the associated 
fluxes of neutrinos for all three bursts. First, the event rates of neutrinos are derived from stan- 
dard fireball phenomenology assuming that the total energy in protons is 10 times more than in 
electrons. In a second approach, we take the existence of the high-energy power-law components 
as indicative of the decay of vr'^-mesons and keep the ratio of the energy going into non-thermal 
electrons and protons as a free parameter. This allows us to calculate the magnitude of the neu- 
trino flux from the related charged pions, and we find that a burst like GRB941017 will produce at 
least one neutrino event and hence be detectable by IceCube if its power-law component extends 
to energies in excess of ^ 50 GeV. While the event rates for the Fermi GST bursts arc small, this 
is due to their large boost factors and rcdshifts as opposed to a burst like GRB941017, which is 
nearby with low boost factor. IceCube and the future observatory KmSNeT will be able to help 
determine if the high-energy components in GRBs are indeed due to photohadronic interactions or 
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Fig. 6.— GRB941017: beamed energy Vs. Cp/ee 



if other scenarios are more viable. 
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